@-PLOS | SUBMISSION 


2105.11503v2 [physics.bio-ph] 26 May 2021 


The Bank Of Swimming Organisms at the Micron Scale 
(BOSO-Micro) 


Marcos F. Velho Rodrigues!, Maciej Lisicki?, Eric Lauga!” 


1 Department of Applied Mathematics and Theoretical Physics, University of 
Cambridge, Cambridge CB3 0WA, United Kingdom. 
2 Faculty of Physics, University of Warsaw, Warsaw, Poland. 


*Email: e.lauga@damtp.cam.ac.uk 


Abstract 


Unicellular microscopic organisms living in aqueous environments outnumber all other 
creatures on Earth. A large proportion of them are able to self-propel in fluids with a 
vast diversity of swimming gaits and motility patterns. In this paper we present a 
biophysical survey of the available experimental data produced to date on the 
characteristics of motile behaviour in unicellular microswimmers. We assemble from the 
available literature empirical data on the motility of four broad categories of organisms: 
bacteria (and archaea), flagellated eukaryotes, spermatozoa and ciliates. Whenever 
possible, we gather the following biological, morphological, kinematic and dynamical 
parameters: species, geometry and size of the organisms, swimming speeds, actuation 
frequencies, actuation amplitudes, number of flagella and properties of the surrounding 
fluid. We then organise the data using the established fluid mechanics principles for 
propulsion at low Reynolds number. Specifically, we use theoretical biophysical models 
for the locomotion of cells within the same taxonomic groups of organisms as a means of 
rationalising the raw material we have assembled, while demonstrating the variability 
for organisms of different species within the same group. The material gathered in our 
work is an attempt to summarise the available experimental data in the field, providing 
a convenient and practical reference point for future studies. 
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A The database of swimming microorganisms 


1 Introduction 


Swimming microorganisms were first observed almost 350 years ago by Antonie van 
Leeuwenhoek [1]. Since then, extensive knowledge has been obtained on their form, 
function, genetics and behaviour EI We now also understand the vital role they play in 
ecosystems |3| as well as in the individual organisms they can inhabit, and whose health 
they influence H Their ubiquity demonstrates an astonishing diversity and 
adaptability to the most extreme conditions. Furthermore, the involvement of swimming 
microorganisms in biological processes, irrespective of habitat, is invariably and directly 
linked to their motility. The chance of a ciliate escaping a predator BE the capacity of 
a spermatozoon to enter and fertilise an egg H and the virulent spreading of pathogenic 
bacteria |8| are but a few examples of how cell motility can be decisive for survival. 

Swimming in a fluid on small, cellular length scales is subject to the physical 
constraints imposed by the viscosity of the fluid. With typical lengths of the order of 
microns, and speeds of a few to hundreds of microns per second, the fluid flows set up 
by microswimmers are characterised by negligibly small Reynolds numbers. The world 
in which their locomotion takes place is therefore dominated by viscous friction and the 
effects of inertia are unimportant [oj]. As a result, the propulsion strategies employed 
by larger organisms such as fish, mammals, insects and birds are ineffective on cellular 
length and time scales [12] [18]. 

Swimming microorganisms have thus developed physical mechanisms to successfully 
overcome, and in fact exploit, viscous drag by actuating slender tail-like appendages 
called flagella [19]. Somewhat confusingly, the same name is used to refer to either the 
polymeric filaments of prokaryotes or the more complex, muscle-like flexible organelles 
of eukaryotes. In the former case, the filaments are semi-rigid and helical, and they are 
rotated passively by molecular motors embedded in the cell wall [29]. For the latter, the 
flagella undergo three-dimensional active motion resulting from the action of 
internally-distributed motor proteins pl. Despite the variation in structure, distribution 
and beating pattern of flagella between species, the actuation of flagella in a viscous 
fluid provides the unifying biophysical picture through which the locomotion of all 
microorganisms can be understood. 

Assessing how fast a certain microorganism can swim is not a simple task. Motility 
is strongly dependent on temperature and on the viscosity of the medium in 
which the cells swim [24428]. Absolute pressure [29], pH and even magnetic 
field have also been shown to influence the motility of certain species. The motile 
behaviour of microorganisms may also change depending on whether they are 
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undertaking the role of prey or predator BEES ES) Furthermore, cellular propulsion 
also depends on biochemical factors . Swimming speeds for different species 
within the same genus (e.g. Vibrio, Ceratium, Peridinium and Paramecium) and even 
different strains of the same species (e.g. Escherichia coli (36][37], Campylobacter 
jejuni and Pseudomonas aeruginosa [38]) are available in the literature but little 
information is given on the variability of the swimming speed within a species or even 
for an individual organism. Overall, data on the swimming speed variability of different 
organisms are rather scarce. Our recent study for eukaryotic microswimmers has shown 
that some of the swimming speed distributions have a universal character when 
appropriately re-scaled but the lack of data limits a more detailed analysis. Since 
motility may be the key factor distinguishing between the regimes of cell feeding 

(i.e. advective vs diffusive) or sensing (e.g. spatial vs temporal) [40], extensive data on 
swimming might aid elucidating the physical mechanisms affecting the cell behaviour. 

The biophysical description of cellular propulsion was pioneered in the last century 
with the works of Gray (from the biology side) and Taylor (mathematics) [42]; and 
it has now grown into a mature field of research HEIEREN Despite many 
theoretical advances, the difficulties of observation and measurement on small scales, as 
well as the complexity of the fluctuating fluid flows continue to offer outstanding 
challenges for detailed studies. In addition, the locomotion of cells links to the rapidly 
growing field of artificial active matter, addressing the question of how microbiology, 
medicine and robotics could work together for the creation and manipulation of 
artificial swimmers, some of which are inspired by flagellated organisms [49]. These 
laboratory swimmers have a promising potential to perform site-specific drug deliveries, 
or chemical sensing, and to assist micro-manipulations in advanced surgery, enhancing 
the effectiveness of medical treatments BOES) 

Motivated by the combination of current activity in the research field and its rich 
scientific history, we carry out in this paper a biophysical survey of the available 
experimental data produced to date (13 April 2021) on the characteristics of motile 
behaviour in unicellular microswimmers. Specifically, we assemble from the available 
published literature empirical data on the motility of four broad categories of organisms, 
namely bacteria (and archaea), flagellated eukaryotes, spermatozoa and ciliates. 
Whenever possible, we gather a broad set of parameters related to biological, 
morphological, kinematic and dynamical aspects of the swimming cells: species, 
geometry and size of the organisms, swimming speeds, actuation frequencies and 
amplitudes, number of flagella and properties of the surrounding fluid. We assemble our 
results in a large downloadable database that we call BOSO-Micro, with BOSO standing 
for “Bank Of Swimming Organisms” and “Micro” emphasising their microscopic scale. 

We then analyse the data from the database in light of the established fluid 
mechanics principles for propulsion at low Reynolds number in order to sort and 
organise the assembled raw material. We reproduce classical scalings for the locomotion 
of cells within the same taxonomic groups, while demonstrating the variability between 
different species within the same group. The resulting database, which is made available 
with this paper and downloadable from the Center for Open Science (OSF) repository, 
provides a convenient and practical reference point for future studies [54]. Despite our 
best efforts, some species and studies may have been left out of our dataset, and since 
research in the field is active and ongoing, it is important to also allow our database to 
be easily and continuously extended. To allow future collaborative effort of the 
community, we have also organised an open source version of the database on 
GitHub [55], which can be supplemented with new data while retaining a version 
control. 

The paper is organised as follows. In Section [2] we describe in detail the structure of 
the database, its sources, and the procedures used for data selection, extraction and 
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Species B U w £ Re Rew 
[um] [ums] [rads '] [um] 

E. coli (bacteria) 25 24.1 823.1 8.3 | 6.7510-5 6.351077 

H. salinarum (archaea) 2.6 33 144.5 4.3 | 9.611076 2.991073 

G. lamblia (flag. eukaryote) | 11.3 26 81.7 11.6 | 3.2810-4 1.22107? 

Bull spermatozoon (Metazoa) | 8.9 97 129.2 54.0 | 9.641074 4.221071 

P. caudatum (ciliate) 242 1476.5 1973 12 | 4.0010-! 3.18107? 


Table 1. Steady (Re) and oscillatory (Re,,) Reynolds numbers for five representative 
organisms from the database. The values of the mass density (p) and dynamic viscosity 
(7) used correspond to water at 25°C. 


processing. We also briefly outline the theoretical basis of locomotion at low Reynolds 
number that serves as a guide for the exploration of our data. We then present and 
discuss the collected data, separating them according to the different taxonomic groups: 
bacteria and archaea (Sec. B). flagellated eukaryotes (Sec. H. spermatozoa (Sec. [5) and 
ciliates (Sec. (6). We summarise the findings in Sec. |7| where we also comment on the 
potential caveats and limitations of our work. We conclude the paper by displaying the 
complete database in Appendix [A] 


2 Methods 


2.1 Propulsion at low Reynolds number 


Cellular swimming is invariably coupled to the fluid mechanics of the surrounding 
environment. Biological locomotion in aqueous media happens on a wide range of 
spatial scales, from sub-micrometre bacteria to whales measuring tens of metres. In all 
cases, steady swimming results from balancing the propulsive forces generated by the 
moving swimmer with the frictional (drag) forces from the surrounding 

environment lot. Propulsion results from the biological actuation, which always 
involves motion of the body relative to the fluid. This in turn generates flow, which 
dissipates energy and thus resists the motion. 

For biological locomotion in Newtonian fluids, the fluid flow around a swimming 
organism is governed by the Navier-Stokes equations. However, in the regime of interest 
for this work, the effects of viscosity on the motion typically dominate inertial effects, as 
classically quantified by the dimensionless Reynolds number. Assuming U to be the 
typical speed scale of a swimmer of a characteristic size B, moving through a fluid of 
mass density p and dynamic viscosity 7, the ratio of inertial to viscous forces is defined 
as the (steady) Reynolds number, Re = pU B/n. Because the propulsion mechanism 
often involves the periodic motion of biological organelles of characteristic length £ and 
angular frequency w, another dimensionless number can be constructed, termed the 
oscillatory Reynolds number and defined as Re, = pwl? /n. 

In Table [I] we estimate both values of Re and Re, for a number of representative 
organisms from the database assuming their environment to be water at 25?C. In the 
majority of cases, these estimates suggest that it is appropriate to neglect all inertial 
effects when compared to viscous forces, as both Re < 1 and Re, < 1, or at most just 
below one. To interpret the dynamics of microswimmers, it is thus appropriate to 
consider the over-damped limit, when the fluid dynamics are governed by the steady 
Stokes equations. For a detailed overview of the fluid dynamics of locomotion at low 


Reynolds we refer to classical work in Refs. 58|. 
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2.2 Data collection and processing 


In this paper we focus on unicellular microorganisms that can swim on their own, either 
using the actuation of flagella and cilia or by periodic deformations of their cell bodies, 
so that they generate net displacements via interactions with the surrounding fluid. We 
therefore do not include gliding and twitching motility, nor amoeboid displacement. 

Swarming bacteria were however included, because swarmer cells are also swimmer cells. 

In order to identify in the available literature the swimming characteristics of 
multiple organisms, we selected six seminal biophysical papers in the field of biological 
fluid dynamics of microscale locomotion (ordered by year of publication): (i) an early 
analysis of microscale swimming by Taylor [42]; (ii) the work of Gray and Hancock on 
the swimming of spermatozoa [59]; (iii) the lecture on the theory of flagellar 
hydrodynamics by Lighthill [56]; (iv) the introduction to life at low Reynolds number by 
Purcell (9); (v) the classical review paper on locomotion by cilia and flagella by Brennen 
and Winet [19]; and (vi) the study on bacterial locomotion in viscous environments by 
Berg and Turner [60]. These papers are commonly viewed by the community as 
groundbreaking biophysical contributions to the field of microswimmer hydrodynamics, 
which is reflected in the number of citations of these works, summing up to over 530] 

In order to construct the database, we first used the Web of Knowledge database to 
assemble two lists of published references: (a) papers that are cited by any of the six 
source papers, (b) papers that cite any of the six source papers. Each of the resulting 
references was then examined to determine whether it contained any measurements or 
reports on the swimming characteristics of any unicellular microswimmer, or if it led to 
other useful references. We acknowledge that our selection of six initial papers is clearly 
biased towards the fluid mechanics and biophysical aspects, yet we hope that by a 
thorough query of the cited and citing papers we managed to sufficiently extend the 
scope of the search to construct a comprehensive and relevant dataset. In order to allow 
further extension of the database to include new and possibly omitted studies, we refer 
to the open GitHub version of it [55]. Note that we reproduce all the collected 
information in the form of tables in Appendix [A] in which we list all relevant material 
in a concise form. 

In addition to the cell swimming speed, we extracted other geometrical and 
kinematic characteristics of the organisms when available in experimental studies. 
'These parameters are summarised on the sketches in Fig. [1] for cells with a small 
number of flagella (top) and for cells with many appendages (bottom): dimensions of 
cell bodies, swimming speeds, lengths and beat frequencies of cilia and flagella, 
wavelengths, wave speeds, amplitudes and form of the propagated waves (two or 
three-dimensional, sinusoidal, helicoidal or complex patterns for flagella, and 
metachrony for cilia [61]). Note that several works exist that review solely the 
morphological features of swimming microorganisms [62][64]. As the focus of our paper 
is on the relationship between geometry, kinematics and locomotion, we chose not to 
include in our database any study that does not report any swimming speeds. 

In all, the database contains a total of 382 species for which we were able to find at 
least one measurement on swimming speed along with other characteristics. Within the 
tree of life, microswimmers of these species are present in all domains: Bacteria and 
Archaea (together encompassing prokaryotic organisms), and Eukaryota (including 
flagellated and ciliated cells and the spermatozoa of multicellular organisms). Members 
of these different groups clearly differ in size, propulsion modes and other physical 
characteristics. In particular, we plot in Fig. [2] the number of flagella (or cilia) of each 
organism against the typical cell body length, demonstrating the partial clustering of 
organisms within their taxonomic groups. On top of variability within taxa, there is a 


1The respective numbers of citations are: 614 [19]; 240 [60]; 733 [59]; 541 [56]; 2461 [o]; 736 [49]. 


Source: Web of Knowledge, 13 April 2021. 
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Figure 1. Top: Geometrical and kinematic parameters of flagellated swimmers, 
illustrated here for a bacterium; we use the same symbols for cells employing planar or 
helical waves for simplicity. Bottom: Geometrical and kinematic parameters of ciliated 
swimmers. Drawings by Marcos F. Velho Rodrigues. 


considerable diversity even within groups, and both parameters can span several orders 
of magnitude. Bearing this in mind, we analyse each taxonomic group separately in 
what follows. 

In order to help visualise the range of the present study, we also follow taxonomy as 
presented in the Open Tree of Life and sketch in Fig. [3] the various phylogenetic 
branches included in our work together with a drawing of one representative organism 
within each phylum covered. 
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Figure 2. Number of appendages, i.e. cilia or flagella, of each organism (whenever 
available) plotted against the cell body length. Both characteristics span orders of 
magnitude but the data cluster within taxonomic groups. 
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Figure 4. Histograms of body lengths, B (um, left), and swimming speeds, U (ums !, 


right), for rod-shaped bacteria (excluding spirochaetes and Spiroplasma) 

((B) = 5.79 + 9.33 um (n = 66), (U) = 48.33 + 98.47 ums"! (n = 77)), spirochaetes 
((B) = 18.59 + 13.02 um (n = 17), (U) = 17.94 + 18.84 ums! (n = 15)), Spiroplasma 
((B) = 5.72 + 0.28 um (n = 2), (U) = 1.69 + 0.81 ums! (n = 2)) and archaea 

((B) = 2.71- 2.12 um (n = 10), (U) = 89.18 + 126.57 um s-! (n = 10)) from our 
database. Most organisms have sizes below 10 um ((B) = 7.75 = 10.85 um (n = 95)) 
and swimming speeds below 100 wms~! ((U) = 46.98 + 95.42 ums! (n = 104)). 


3 Bacteria and Archaea 


We start our journey through swimming microorganisms with prokaryotes, namely the 
domains Bacteria and Archaea. Bacteria constitute the bulk of the biomass on Earth, 
inhabiting the soil, water reservoirs, and the guts of larger organisms. They are simple 
cells without a nucleus, yet they display a remarkable diversity of shapes [66]. Motility 
is a crucial feature for many species of bacteria, in particular for nutrition purposes, and 
to this end bacteria have developed various propulsion strategies [67]. 

Two broad categories of swimming bacteria exist. In the first one, propulsion is 
enabled by the actuated motion of flagella located in the fluid outside the cell body [20]. 
Unlike their active eukaryotic analogues, prokaryotic flagellar filaments are passive 
organelles of typical length of a few microns, attached to a flexible hook that acts 
as a joint connected to a molecular motor embedded in the cell wall. The word 
flagellum (plural flagella) is used to refer to the motor-hook-filament complex. The 
bacterial rotary motor, driven internally by ion fluxes, exerts a torque on the hook, 
which transmits it to the filament thereby inducing its rotational motion. Because the 
flagellar filaments have helical shapes, their rotation in a viscous fluid induces a 
hydrodynamic propulsive force and leads to the motion of the organism [10]. 

Flagellated bacteria can be equipped with anything from one flagellum 
(monotrichous cells) to a few flagella originating from different points on the cell body 
[69]. Polar bacteria have their flagella positioned in the vicinity of the pole of the cell. 
Other arrangements are seen in lophotrichous (a tuft of flagella at the pole) and 
amphitrichous (flagella at each pole) cells, while for peritrichous species (including the 
well-studied model organism Escherichia coli) the rotary motors are located 
approximately randomly on the cell body. 

Some species of flagellated bacteria can also display a mode of motility named 
swarming, where cells undergo changes in morphology and rely on intercellular 
interactions to move near surfaces [79]. Some species can transition from swimming to 
swarming behaviours by relying on polar flagella for swimming, while exploiting several 
flagella distributed along the sides of their bodies for swarming [71]. 'The data for most 
bacteria in our database is presented in Table [4] 

In the second type of bacterial swimming, cells move via a time-dependent 
deformation of their body. Famously, cells in the phylum Spirochaetes are 
morphologically distinguished by having internal axial flagellar filaments running 
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Figure 5. Histograms of aspect ratios W/B (left) and body-to-flagellum length B/L 
(right) for rod-shaped bacteria (excluding spirochaetes and Spiroplasma) 

((W/B) = 0.33 + 0.20 (n = 63), (B/L) = 0.93 + 1.19 (n = 28)), spirochaetes 

((W/ B) = 0.02 + 0.01 (n = 17)), Spiroplasma ((W/B) = 0.03 + 0.00 (n = 2)) and 
archaea ((W/B) = 0.11 + 0.06 (n = 2), (B/ L) = 0.63 + 0.24 (n = 9)). All bacteria in 
our study are prolate, with an average aspect ratio (W/ B) = 0.25 + 0.22 (n = 84), with 
a notable slenderness of spirochaetes and Spiroplasma. If the prokaryotes possess freely 
rotating flagella, their length often exceeds the body size (B/L) = 0.86 + 1.05 (n = 37) 
(both spirochaetes and Spiroplasma are not included in the B/L graph). 


lengthwise between the inner and outer membrane of their periplasmatic space, 
producing helical waves in the cell body with no apparent slippage with respect to the 
surrounding fluid [72]. Unlike typical rod-shaped bacteria, this particular configuration 
allows them to swim in extremely viscous gel-like media. 

Finally, cells in the genus Spiroplasma do not present axial flagellar filaments. 
Instead, they swim by propagating kink pairs along their helical body using the motion 
of its cytoskeleton. This creates a processive change in the helicity of the body, which 
also allows them to move through extremely viscous fluids [73]. Our data for 
spirochaetes and Spiroplasma is presented in Table [5] 

Relatively less studied are the species in the prokaryotic domain Archaea. Archaea 
have morphologies similar to bacteria but, equipped with a different molecular 
organisation, they are able to live under conditions that are extreme and hostile to 
other forms of life. Other differences exist; for example, some species of archaea have 
square-shaped bodies, unlike any bacterium or eukaryote Glo Although the 
actuation of archaeal flagella has been characterised in detail , the motile behaviour 
of only about 10 species in the whole domain has been studied so far, with all data 
summarised in Table [6] 


3.1 Geometry and swimming speeds of the cells 


The distribution of sizes and speeds of the prokaryotes from Tables [4] [5] and [6] are shown 
in Fig. 4| The characteristic length of the cell bodies does not exceed 10 wm while the 
typical swimming speeds are of the order of tens of ums !. 

'The shapes of the prokaryotes are next quantified in the distributions shown in 
Fig.|5| (left). The cells are close to ellipsoidal, with an aspect ratio W/B (body width to 
length) not exceeding 1 and an average of about 0.25. In contrast, spirochaetes and 
Spiroplasma are slender, with the aspect ratio not exceeding 0.05. We also plot in Fig. B] 
(right) the distribution of body-to-flagellum lengths for cells with external flagellar 
filaments (i.e. excluding spirochaetes and Spiroplasma). This is typically smaller than 
unity, indicating that the helical filaments are longer than the cell body in most cases. 

'The swimming speed for all prokaryotes in our database is plotted in Fig. [6] against 
the cell body length (top panel) and width (bottom panel), with colours used to divide 
our dataset into the specific taxonomic groups. Clearly, a wide spread of values exist for 
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Figure 6. Swimming speed, U (ums~'), as function of the cell body length, B (um, 
top), and body width, W (um, bottom), for all our registered prokaryotes. Error bars 
represent standard deviations, whenever available, or the span between the recorded 


maximum and minimum values. 


the swimming speeds and in the next section we use a mathematical model for bacterial 
locomotion in order to gain further insight into the data. 


3.2 Modelling of swimming for flagellated prokaryotes 


We focus in what follows on the case of rod-shaped prokaryotes. The flagellar 


locomotion of bacteria relies on the motor rotation being transmitted to the passive 
flagellar filament via the flexible hook [20]. 'The rotation of the motor is generated by 
ion fluxes and in the forward propulsion mode the rotary motor works at approximately 
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constant torque [77]. The value of this torque, however, has been under some debate. 
Berry and Berg estimated the stall torque in an optical tweezers experiment to be of the 
order of 4600 pN nm [78], while Reid e£ al. attached micrometer beads to flagella to 
measure the motor torque to be 1260 + 190 pN nm [79]. In magnetic tweezers 
experiments involving the attachment of paramagnetic beads, the corresponding torque 
amounted to 874 + 206 pN nm [80]. In contrast, a simplified theoretical model predicts 
a lower value of 370 + 100 pN nm while recent numerical simulations reported 
values in the range 440 — 820 pN nm [82]. Kinosita et al. managed to observe in 
detail the flagellar rotation of the archaeon Halobacterium salinarum and estimated its 
motor torque to be as low as 50 pN nm. However, different species of bacteria can have 
very different motor structures [83], which leads to a wide range of possible values for 
the propulsive torque [84]. 

In order to estimate the motor torque of various species in our dataset, we consider a 
generalised mathematical model for the swimming of flagellated prokaryotes. For 
simplicity we focus on the case of a cell rotating a single helical filament [85]. However, 
the resulting model should remain valid for a prokaryote with many flagella, since 
during swimming all flagellar filaments bundle together and rotate as if they formed a 
single helix [11]. Furthermore, as we show later, the model can be easily adapted to 
cope with the impact of bundled flagella. 

A prokaryotic flagellar filament of length L is well approximated by a rigid helix of 
pitch A and radius h (as shown in Fig.|1| top), rotating at an angular velocity w = 2r f 
relative to the cell body, where f is the frequency of rotation of the flagellum. In turn, 
the cell body rotates at an angular velocity €) relative to the fluid to maintain an overall 
torque balance on the cell. At low Reynolds number, the helical filament is subject to a 
hydrodynamic thrust F and a viscous torque T which depend linearly with the axial 
speed U and the rotation rate of the filament relatively to the fluid Q +w. This may be 


written as 
flagellum f f22 W i 


Prokaryotic flagellar filaments are very thin, with typical radii of 0.02 wm and 
average lengths a thousand times greater (the mean value of all lengths in our database 
is (L) = 8.28 um), so that the coefficients of the symmetric matrix fi; can be evaluated 
using the resistive-force theory of viscous hydrodynamics valid for very slender filaments 
[59]. Integrating the local hydrodynamic forces on each small segment of the flagellum 
using the viscous drag coefficients per unit length, c} and cj, in the directions normal 
and tangential to each segment respectively (see details below), yields the classical 
result that the resistance coefficients in Eq. are given by 


fu = (e cos? 0 + c, sin? 0) L, (2a) 
fiz = (ci-—cj)sin0cos8 hL, (2b) 
fa = (ecos? + Cj sin? 0) h2 L. (2c) 


where 0 = arctan(27h/A) is the helix pitch angle DES) 

The cell body, modelled as a prolate spheroid of length B and diameter W (Fig. 
top), is subject to a hydrodynamic force F' proportional to the swimming speed U and a 
hydrodynamic torque T proportional to its rotation rate Q. Assuming the cell to rotate 
about its principal axis leads to 


(r). -- (o J (3) 


where the off-diagonal coefficients vanish due to the symmetry of the body. 
During steady, straight swimming, the sum of forces and torques on the swimming 
bacterium must vanish, and thus combining Eqs. and we obtain a linear system 
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of equations for the swimming speed and angular rotation as a function of the rotation 
rate of the filament as 


bir + fu fia U fiz 
= — . 4 
( fia b22 + f22 Q f22 * E 
Solving for U and Q as functions of w leads to the relations 
U = fizb22 w (5a) 


f — (bii + fi) (b22 + bal ` 
= foo(fir + 611) — fÓ " 
Boe ac c D Sd) ^ 2 


The torque Tm exerted by the flagellar motor is, by definition, given by 
Tm = fi2U + Gatti + w), which after substitution into Eq. yields 


boa (feo — f22(b11 + fui)) 


Tu = w, 6 
™ fi — (bir + fir) (b22 + f22) (6) 
and therefore the ratio between the swimming speed and the torque exerted by the 
motor is only a function of the various resistance coefficients, as 
U _ fiz (7) 
Tm f — fə2(bui + fu) 


The ratio between f2, and fiif» can be computed using the expressions given by 
Eq. and we obtain 


fi (c1 — eq)? sin? 0 cos? 0 


= (8) 


fufe (cl cos20 + ci sin? 0) (cy sin? 0 + cÀ cos? 6) 


The right hand side of Eq. is always positive (since c) cj > 0). Its derivative with 
respect to 0 is given by 


2c, ole) — cy)? sin 0 cos 0(cos2 0 — sin? 0) 


(c1 sin? 0 + eq cos? 0)2(ce | cos? 0 + cy sin? 0)2 ' 


(9) 


which has 0 = (kz /4,k € Z} as roots, for all values of cj and c1. Since 0 = (0,7/2) are 
zeros of Eq. themselves, 0 = 7/4 gives the maximum possible value for the ratio as 


fà -la-a Gast, 


fir fo = (cr +o)? Xl ae 


(10) 


It is usually a good approximation to take cy/c, 7 1/2, so that the ratio f?5/(fi1f22) is 
bound from above by 1/9, and one may thus approximately neglect the contribution of 
f£; in the denominator of Eq. (7), yielding the simplified result 


o fe 
foo(bir + fu) 


The drag coefficient b11 for a prolate spheroid of length B and diameter W depends 
on a geometric factor Crg that involves the eccentricity e of the spheroid, given by 


e= VI (W/B)2 (0 < e <1), as [87] 


U 
Tm 


(11) 


(12) 


8 1 T 
bi; = 307BCra(W/B), CrB= ~ Län (1 + e?) log u 


1—e 
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Figure 7. Propulsion speed of rod-shaped prokaryotes vs morphological factor 1/£?. 
Bacteria are plotted in squares and archaea in circles with colours used to distinguish 
between the different taxonomic classes. The plot, along with Eq. (14), allows to 
estimate the range of bacterial motor torques 27.48 — 1907 pN nm, represented by the 
shaded area. 


'The asymptotic limit of very slender spheroids, evaluated in Ref. , also gives the 
friction coefficients for the motion of a rod of length L and maximal thickness 2b as 


Jon 270 


^L 7 log(L/b) +1/2' “| log(L/b) — 1/2’ 


(b/ L < 1), (13) 
which, for large aspect ratios, yield the approximation cj /c, œ% 1/2, as above. Assuming 
for simplicity the pitch angle to be 0 ~ 7/4, and using the friction coefficients as in 


Eq. (13), Eq. takes the final explicit form 


Tm 
=: 
where the characteristic length £ depends solely on the morphology of the swimmer and 
results from the interplay of body and flagellum size. The result in Eq. relates 
therefore the swimming speed U to the flagellar motor torque Tm via the viscosity of 
the fluid (7) and a morphological factor (€). Note that by adjusting the helix thickness 
2b, the model can address the impact of having several filaments inside the flagellar 
bundle [58]. Since the effect of b in Eq. is logarithmic, its impact on our results is 


minimal. 


U € = y9rh (BCrs (W/B) + [L/ (log(L/b) + 1/2)]), (14) 


3.3 Insights from data 


We can now use the model introduced above in order to help organise our database and 
provide a simple estimate of the range of motor torques in the available data. In Fig. 
we plot the swimming speed, U, for rod-shaped bacteria and archaea as a function of 

the morphological factor 1/£? for all the species for which our database gives access to 
the parameters involved in the definition of € in Eq. (we assumed the thickness of 
the flagella to be 2b = 0.02 um in all cases). The ratio between U and 1/£? should yield 
an estimate of the effective flagellar motor torque, Tm. An important limitation is that 
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the value of the viscosity is, alas, rarely given directly in the studies gathered in our 
database. We thus assume the viscosity 7 in Eq. to be that of water at 25°C and in 
Fig. [7| we display the range of motor torques so obtained using parallel lines enclosing 
the shaded area. The lower and upper bounds of the motor torque Tm are obtained to 
be 27.48 pN nm (for Halobacterium salinarum) and 1907 pN nm (Pseudomonas 
fluorescens). This large range highlights the intrinsic variability within this group, 
corresponding to the observed scatter of the data. 


4 Flagellated eukaryotes (excluding spermatozoa 
and ciliates) 


Eukaryotic cells are not just morphologically distinct from prokaryotes, they also have 
different important biological features, including the presence of a cellular nucleus. 
Their propulsion is enabled by an internal mechanism that is fundamentally different 
from, and more complex than, that of prokaryotes. The central structure of eukaryotic 
flagella and cilia is termed the axoneme and is usually composed of nine microtubule 
doublet filaments surrounding a tenth central pair of microtubules. Cross-linking dynein 
motors allow the relative sliding of the microtubules, which results in the propagation of 
bending deformations along the flexible flagellum that can take the form of 
travelling waves, either planar or helical, as well as of complex two- (2D) and 
three-dimensional (3D) kinematics. 

The eukaryotic flagellar waves usually propagate from the flagellum base to its distal 
end, but some species have waves travelling in the other direction. Similarly, while most 
species swim with flagella trailing, some, such as the alga Ochromonas danica, 
self-propel with their flagella leading the cell. We refer to Jahn & Votta for an extensive 
overview of the observed beating patterns [90]. One of the most fundamental beating 
patterns displayed by eukaryotic cells is a simple planar sine wave, and we will use it as 
the basis for the modelling introduced below. Note that flagella of some eukaryotic 
species display perpendicularly attached rigid structures, termed mastigonemes, which 
give a hairy microstructure to the flagellum and allow the cells to generate propulsion in 
the same direction as the propagating wave [19]|p1]. Some algae such as 
Chlamydomonas do not even rely on waves to swim, but do so by swinging a pair of 
short flagella in a breaststroke way. 

Eukaryotic cells are generally one or two orders of magnitude larger in size than 
prokaryotes and are therefore more easily observed experimentally. A number of past 
review papers gathered swimming speeds and body lengths for tens of 
organisms [92][95]. Our database builds on, and extends, these datasets by introducing 
a number of new important cellular parameters and new organisms. Note that parts of 
our data for eukaryotic cells, particularly the average sizes and swimming speeds have 
been published elsewhere [39]. 

Among swimming unicellular eukaryotes, three families with different morphology 
can be distinguished: flagellates, spermatozoa, and ciliates. Flagellates — the focus of 
this section — typically possess a few long flagella attached to their bodies, which they 
actuate in order to achieve propulsion (for organisms in this section, the typical number 
of flagella rarely exceeds 10). Spermatozoa are also remarkable flagellated swimmers but 
they lack the ability to reproduce, thus are not considered living organisms. Lastly, 
ciliates are much larger in size and are covered by dense arrays of cilia, which are short 
flagella that move collectively to create flow along the cell surfaces. The qualitative 
difference in their swimming speeds, as well as their geometric characteristics such as 
their size and their number of flagella, warrants separate statistical analysis for each 
group [39]: spermatozoa are therefore addressed in Sec. [5] and ciliates in Sec. [6] 
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Figure 8. Histograms of body lengths, B (um, left), and swimming speeds, U (ums !, 


right), for flagellated eukaryotic swimmers (excluding spermatozoa and ciliates) in our 
dataset. The average cell length is (B) = 38.87 + 56.64 um (n = 113) and the average 
swimming speed (U) = 186.70 + 208.77 ums^! (n = 116). 
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Figure 9. Histograms of aspect ratios W/B (left) and body-to-flagellum length ratios 
B/L (right) for flagellated eukaryotic swimmers. For all organisms in this category, the 
aspect ratios do not exceed ~ 1.1, and the shape distribution indicates a slightly prolate 
shape on average, with (W/B) = 0.60 + 0.27 (n = 73). The distribution of 
body-to-flagellum length ratios shows that flagella tend to be of length comparable to 
the cell body, with a few exceptions (B/L) = 1.03 + 0.79 (n = 49). 


4.1 Geometry and swimming speeds of the cells 


'The typical sizes and swimming speeds of eukaryotic flagellates are presented in Fig. 
based on the data from Table [7] Significantly larger and faster than prokaryotic cells, 
the distributions are dominated by the low-values tails. 

Most cells are close to the average values, with several outliers in the large size and 
speed ranges. The statistical properties of these distributions have been discussed in 
detail in our previous work [39]. We may gain further insight by considering the 
distribution of aspect ratios for the cell bodies, W/ B, and the relative cell 
body-to-flagella lengths, B/L, both of which are shown in Fig. oa The wide distribution 
of aspect ratios confirms that most flagellates are slightly prolate, although several more 
elongated swimmers are also reported. In addition, for most cells the ratio of body to 
flagella length does not exceed 1, confirming that the length of the flagella is 
comparable to the cell size. This feature allows to distinguish flagellated eukaryotes 
from spermatozoa and ciliates. 

In Figs. [10] and [11] we next show how the swimming speeds U of the flagellated 
eukaryotes in our database vary with the flagellar beat frequencies f and flagellar 
lengths L, respectively. Both plots show large variations and no clear trend is evident. 
In the next section we will then adapt the classical derivation by Gray & Hancock 
as a minimal model for the propulsion of eukaryotic flagellates to see the role played by 
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Figure 10. Swimming speed, U (ums), plotted versus the frequency of flagellar 
beat, f (s~'), for flagellated eukaryotes in our dataset (excluding spermatozoa and 
ciliates). Colours mark different classes and sub-classes. Wave-producing organisms are 
plotted in squares and the remaining flagellated eukaryotes are plotted in circles. 


these (and other) parameters in eukaryotic propulsion. 


4.2 Modelling of swimming for flagellated eukaryotes 


We base the description of the locomotion of flagellated eukaryotes on the assumption 
that swimming results from planar travelling waves induced in one or more flagella, 
which push a spheroidal cell body forward. 

The shape of the wave is described in Cartesian coordinates by y = y(z,t), where x 
is the direction of cell movement (see Fig. [12). An infinitesimal segment of the flagellum 
of length ós inclined at an angle 0 to the axis of movement e, is then subjected to a 
hydrodynamic force perpendicular to its orientation, and given by 


ôF, = cı (U, cos 0 — U sin 0) ôs, (15) 


and to a force tangential to the segment given by 
ôF = c (Uy sin 0 + U cos 0) däs, (16) 


Here U and U,(z,t) are the local velocities of the flagellum relative to the fluid in the 
directions along and perpendicular to the overall direction of cell motion, respectively. 
Furthermore, similarly to the section on prokaryotes, c; and cj are the drag coefficients 
per unit length in the directions normal and tangential to ds, respectively (see Eq. [13). 

These two force components produce an infinitesimal net thrust along the x 
direction, ôF = 6F sin 0 — dF) cos0, which we rewrite as 


(cı — c)Uy tan — U (cy + c, tan? 0) 


Ab = 
1+ tan? 0 


As, (17) 


Taking into account the normal speed to be U, = Oy/0t, using tan 0 = Oy/Ox and 
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Figure 11. Swimming speed, U (ums™t), vs length of flagella, L (um), for flagellated 
eukaryotes in our database (excluding spermatozoa and ciliates). Taxonomic classes are 
marked by colours. Wave-producers are again plotted in squares, while other flagellates 
are plotted in circles. 


ôs? = ĝy? + óz2, we transform Eq. into 


T: (18) 


We now need to specify a particular wave form of the beating pattern. One that is often 
observed in eukaryotic swimmers is a planar travelling wave which we approximate 
by a single sine wave of fixed amplitude h, wavelength À and beat frequency f 


2c 


y(z,t) = hsin ( CS d , (19) 


where c = Af is the speed of the propagating flagellar waves. Substituting the sine wave 
into Eq. (18), and taking the slender limit c, ~ 2c), yields 


cA? — U(1 + 242) 
T, 
it A2 


where A = 0y/Ox = (2zh/ À) cos (2n(x + ct)/A). It is convenient to introduce the 
number of complete waves nw in the flagellum of length L, defined as 


1 ye 1 fò 
—=7/ i-r] v 14 A? dz. (21a) 


F = oi (20) 


Because the integrand v1 + A? is a function of period À, a simple substitution shows 
that the number of waves is constant in time, and is given by 


1 À Qh 
= A 21b 
= = (5. (210) 
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Figure 12. Sketch of a swimming eukaryote (spermatozoon of Chaetopterus, Annelida) 
propelled by a single flagellum. We distinguish a section of length ós inclined at an 
angle 0 to the direction of motion e,, which we use to determine the local 
hydrodynamic forces exerting on the flagellum. Drawing by Marcos F. Velho Rodrigues. 


where the auxiliary integral A is 
1 2m 
A(a) = — V 1 + a? cos? a da. (22) 
2m Jo 


With the net thrust ôF in Eq. being also of period À, a good approximation of 
the total thrust produced by the entire flagellum independent of time is given by 


À 
Nw ôF = NCA (en Ea UI» (=) , (23) 
z-—0 


where we have introduced the two auxiliary integrals 


2 
1 2" aq cos2 o 1 T 1+ 2a? cos? a 


= da, Ig(a)= da. 
2T Jo V1+ a? cos? a ala 2T Jo V1+ a? cos? a 
The three functions A, T, and I> are easy to evaluate numerically. Alternatively, by 


writing cos? a = (1 + cos 2a)/2, and neglecting the contributions of the terms in cos 2a 
in the expressions of Eqs. and (24), one gets explicit approximations 


D(a) (24) 


Lia 


ef 
ER h(a) x ——— h(a) = ——— 
E disc 

2 2 


Numerical evaluation of the exact expressions for A, T, and I> shows that the 
approximations above hold to within 13% accuracy for all values h/A < 1. 

For the sake of simplicity, we shall suppose that an organism with N beating flagella 
is subject to a total thrust equal to N times the thrust generated by each flagellum and 
given by Eq. (23). We therefore neglect hydrodynamic interactions between the flagella, 
which we assume all beat collinearly along the swimming direction. 

Steady swimming requires the thrust produced by the flagella to be balanced by the 
drag acting on the cell body. The latter is modelled as a prolate spheroid of length B 
and diameter W. The balance of forces acting on the microorganism along z is then 


(25) 


À 
Nn, | 5F —30nUBCrp = 0, (26) 
p—0 
with Crg(W/B) given by Eq. (12). The swimming speed U can thus be written as 

2th 

j 2 
= ( x ) (27) 

Af ` 1 (28h) , 30nBCre(W/B)’ 
a À l Nnwc| À 
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Figure 13. Swimming speeds of flagellated eukaryotes (excluding spermatozoa and 
ciliates) reported in the database plotted against the theoretical prediction of Eq. (28). 
Colours mark different classes. Square symbols mark organisms for which the prediction 
was directly calculated from the available data, while circles represent organisms for 
which either the body width or one of the flagellar characteristics has been estimated 
(see text for details). 


or using the definition of n,, in Eq. as 


27h 
U ` DU: (28) 
Af ` 2xh 329 B 2nhN ` 


By using o = 277 [log(L/b) — LIT as in Eq. (13), and approximating integrals I4, I2 
and A with the expressions in Eq. (25), we arrive at the final expression 


94252 Aer Zb D 242? L 1 -1 
Hl, E eene 2") (los S . e» 


À X? 2N À b 2 


4.3 Insights from data 


We can now use our model to help organise our data on flagellated eukaryotes. In 

Fig. we compare the swimming speeds from our dataset with those predicted by the 
theoretical model in Eq. (28). Square symbols mark organisms for which all the 
quantities needed to calculate the predicted speed were available. The species plotted in 
circles in the figure had their data incomplete. Whenever the body width W was 
unavailable, we estimated its value using the average aspect ratio (W/B) = 0.60 of 
Fig. [9] When one parameter of the flagellar wave was missing, we estimated it with the 
help of Eq. (21). The radii of the flagella were all fixed at b = 0.2 um. 

In Fig.|13| we see a cluster of data points (mostly the class Kinetoplastea) that 
correlate well with the expected linear dependence. However, many of the organisms 
have a swimming speed that significantly exceeds the predicted values. This may point 
to other mechanisms being involved, such as different beating patterns and cell body 
shapes, which would require a more careful examination. Nevertheless, the basic 
framework proposed by the model provides a useful estimate of the lower bound for the 
swimming speed, which can be exceeded by adopting more effective locomotion 
strategies suited to the organism and its environment. 
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Figure 14. Histograms of body lengths, B (um, left), and swimming speeds, U 
(ums™t, right), for the spermatozoa in the database. The average cell length is 
(B) = 12.21 + 17.25 um (n = 39), while the the average swimming speed is 

(U) = 127.23 + 78.49 ums! (n = 52) over a wide distribution. We use colours to 
distinguish between the different taxonomic classes. 


5 Spermatozoa 


The motile behaviour of the spermatozoa of animals has been studied in detail since the 
beginnings of microscopy due to its importance for reproductive health. Because a 
correlation between motility and fertility has been shown to exist |96/97], numerous 
species of fish [98], birds [99], mammals [A1|[1oo]101], insects and sea 
urchins have had their spermatozoa examined. A particular focus is often placed 
on the relation between either the swimming speed or the amplitude of lateral 
displacement of the cell body and the success in fertilisation by human spermatozoa [7]. 

A remarkable geometrical characteristic of spermatozoa, at least in comparison with 
other flagellated eukaryotes, is their relatively small heads compared to the length of 
their flagella. Despite this difference, the flagella of spermatozoa have the same 
structure detailed above for other eukaryotic cells, and are likewise capable of creating 
complex waveforms. The mathematical modelling of flagellar locomotion outlined in the 
previous section is thus also applicable to the case of swimming spermatozoa. 

Our database of swimming spermatozoa contains 60 different species, for which 
various geometric and dynamic data were found. These include sperm cells of the 
taxonomic classes Insecta, Actinopterygii, Mammalia, Amphibia, Polychaeta, 
Ascidiacea, Echinoidea, Aves, and Bivalvia. As mentionned above, databases of 
morphological measurements for over 400 spermatozoa, particularly of mammalian 
species, are available in literature but since they do not include motility data 
they are not included in our database. 


5.1 Geometry and swimming speeds of the cells 


The distribution of cell body sizes and swimming speeds of spermatozoa are shown in 
Fig. [I4] based on the data from Table[8] With body sizes hardly exceeding 30 m 
(except for one outlier, the cricket spermatozoon, with a size of over 100 jum), we see 
that spermatozoa are typically small compared to other eukaryotic cells. The 
distribution of swimming speeds is relatively uniform, reaching up to 300 uyms~!. While 
their average speeds are close to those of flagellated eukaryotes from Sec. |4| the 
distribution of speeds is dramatically different, deviating from the log-normal seen for 
other flagellated eukaryotes [39]. 

A further inspection of the geometry reveals that the distribution of sperm cell 
aspect ratios (Fig. left) is widely spread, ranging from elongated to spherical. A 
clear distinguishing feature for spermatozoa is the body-to-flagella length ratio (Fig. 
right), which is peaked at small values, showing that the spermatozoa of most species 
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Figure 15. Histograms of aspect ratios W/B (left) and body-to-flagellum length B/L 
(right) for spermatozoa (colours mark the different taxonomic classes). The distribution 
of cell aspect ratios is rather wide, and yields an average value of 

(W/B) = 0.47 + 0.30 (n = 31). The size-to-flagellum length ratios are mostly close to 
the average (B/L) = 0.17 + 0.18 (n = 38), showing that in spermatozoa the flagellum 
length is typically much larger than the cell body. 


= 


have flagella that are over fivefold longer than their body sizes. 


5.2 Hydrodynamic model for locomotion 


The locomotion of flagellated spermatozoa follow the same hydrodynamic principles as 
discussed in detail in Sec. |4| We may thus use as our starting point the the result in 
Eq. (28), which upon using the drag coefficient c = 2z7[log(L/b) — 1/2]! and N = 1 
takes the form 
5 _ = I, (27h/d) f (30) 
I, (2zh/ À) + ap FB OW/ B)AQmh/ X) [log(L./0) — 1/2] 


Note that the second term in the denominator of the right-hand side of Eq. is the 
hydrodynamic load of the dragging cell body, which we include although the flagella are 
notably longer than cell bodies for spermatozoa. 


5.3 Insights from data 


We again turn our attention to the behaviour of the swimming speeds for the cells. In 
Fig. we examine the dependence of the spermatozoa swimming speed U on the 
flagellar beat frequency, f. With most spermatozoa operating in the frequency range 
between 10 and 100 Hz, and swimming speeds of up to 300 uum s^! , we observe a 
pronounced correlation between these two variables across our database. In Fig. we 
also show the dependence of the swimming speed U on the flagellar length L, which 
ranges from about 20 to 120 um. Here, in contrast, no direct or apparent correlation is 
seen between U and L. 

'To help organise the information on the locomotion of sperm cells in our database, 
we resort to the model from Eq. (30), which we compare with the collected data in 
Fig. Circles represent organisms for which either the body width W was unavailable 
(in which case we assumed W/B = 0.47, the average value from Fig. IE or for which 
one parameter of the flagellar wave was missing (and was thus estimated using 
Eq. (21)). The thickness of the flagella was fixed at 2b = 0.4 um. We see that the 
model of Eq. is able to capture the essence of spermatozoan swimming, and better 
than it did for flagellated eukaryotes in the previous section. The outliers can likely be 
explained by the use of more complex wave patterns in some species. 
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Figure 16. Swimming speeds, U (ums-1), as function of flagellar beat frequency f 
(s71), for spermatozoa. A strong correlation between U and f is apparent on the figure. 


6 Ciliates 


Within the diverse group of flagellated eukaryotes, the final family of organisms is 
distinguished by their remarkably large number of flagella, ranging from hundreds to 
tens of thousands (see the distribution in Fig. Ë). These flagella are short compared to 
the size of the cell body and are called in this case cilia — hence the name of ciliates 
given to the whole group. Ciliates have developed a locomotive strategy relying on the 
phased beating of their many cilia. Typically, a single cilium beats using a two-stroke 
motion with a power stroke of an extended cilium followed by a recovery stroke where 
the cilium is curved, generating a polarised beat [10]. From the phased beat of 
neighbouring cilia, collective motion is induced that pumps the surrounding fluid [19], 
thus creating the hydrodynamic forces necessary for locomotion. This collective 
sequential movement of cilia is often observable through the so-called metachronal 
waves of deformation travelling over the surfaces of ciliated cells, resembling spectator 
waves in stadiums. Yet, the underlying ciliary structure is not easily observed and only 
a few studies report successfully the wavelengths of metachronal waves and ciliary beat 
frequencies. In particular, for the model organisms in the genus Paramecium the 
frequencies of ciliary beat of all the different regions of the cell have been accessed [107]. 

The mathematical modelling of metachronal waves can be undertaken at various 
levels of complexity [10/12], starting with coarse-grained continuum models, such as the 
squirmer model , up to detailed simulations of the deformations of individual 
cilia interacting hydrodynamically ol. Non-hydrodynamic interactions via 
intra-cellular coupling mediated by the cell body are also important (112)[113]. 
Independently of the specific coordination mechanism, ciliates all swim by transporting 
the surrounding fluid along their surfaces, and move in the direction opposite to the 
fluid motion. By using different models for this effective transport mechanism, we can 
now test several hypotheses across our database of ciliates, which involves data for 93 
species. Note that the distribution of swimming speeds across species from this dataset 
has been published in our earlier contribution [39]. 


6.1 Geometry and swimming speeds of the cells 


In Fig. we present histograms of sizes and swimming speeds for the ciliates in our 
database. Most of the organisms are close to, or slightly below, average values, which is 
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Figure 17. Swimming speeds, U (ums~'), as function of flagellar lengths, L (um), for 
spermatozoa. In contrast with the result in Fig. no clear correlation between U and 
L is observed here. 


highlighted by the skewness of the distributions [39]. The cells are notably larger 
(average length about 200 um) and faster (average speed of over one millimetre per 
second) than any other group in our database. As a result, the dimensionless Péclet 
number for relevant molecular solutes (such as ions) around the ciliates is of the order of 
100 which means that, in contrast to bacteria and flagellates, ciliates live in a physical 
environment where advection and thus the ability to stir the surrounding fluid may be 
the life-driving mechanism [39]. 

The distribution of aspect ratios of the cells, along with the body-to-cilia lengths, 
are shown in Fig. 'The former peaks at the mean value of about 0.5, indicating 
prolate cell bodies. The large values of the body-to-cilium length ratios confirm that 
cilia take the form of tiny hairs covering the cell body, much smaller than the body 
itself. This in turn justifies coarse-grained modelling approaches representing the cell 
body as a continuous surface capable of exerting stress, thereby locally averaging the 
collective motion of individual cilia. 


6.2 Models for ciliary propulsion 


In search of means to organise our data on the locomotion of ciliates, we propose below 
three distinct ciliary propulsion models that each assume a different property to be 
constant among the cells during forward swimming. These three approaches model the 
swimming of the cells as induced by: (A) a constant tangential stress exerted on the 
cell surface by the cilia array; (B) a constant force exerted by each individual cilium on 
the fluid; (C) a constant effective fluid speed induced near the cell surface by the cilia. 
We model a ciliate cell as a prolate spheroid of length B and diameter W. We set the 
x-axis along the long axis of the cell, taken to also be the direction of movement. The 
ciliate swimming with speed U along x is then subject to a viscous drag of magnitude 


D= 3rnBC FBU, (31) 


with the geometry-dependent coefficient CF g in Eq. (12). Balancing this drag with the 
propulsive force generated by the collective action of the cilia yields different models for 
the swimming speed U, according to how one exactly models the propulsive force. 
Some aspects of the mathematical description of the cell will be useful in what 
follows. A cross-section of the spheroid containing e, is an ellipse of eccentricity 
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Figure 18. Reported propulsion speed of spermatozoa compared with the values 
predicted by the theoretical model in Eq. (30). Colour scheme distinguishes between the 
different taxonomic classes. Squares represent spermatozoa that had all parameters 
available in the literature, while the circles mark cases where at least one parameter had 
to be estimated (via (W/B) = 0.47 from Fig. [15] or through Eq. Ei 


e = /1— (W/ B)?. Every point of the ellipse can be parametrised in polar coordinates 
by 
W 


2/1 — (ecos0)?' 


with the origin placed at the centre between its foci. Every point on the surface of the 
spheroidal body can then be written using spherical coordinates as 


r(0) = (32) 


x(0, p) = r(0) [cos0 e; + sin 0 (cos pe, + singe;)], 0 € [0,7], v € [0, 2r]. (33) 


One may thus write the axisymmetric, unit vector tangential to the spheroidal surface 
and pointing along the polar angle as 


l '(0)e,. + r(0)e 
t(0) = VOZIO [r (0)e, + r(0)eo] , (34) 


where r'(0) = dr(@)/d@. Finally, an infinitesimal surface element on the spheroidal 


surface is given by 
dS = r(0) sinf /r(0)? + r (0)? dé dy. (35) 


Let then x be a given point on the spheroidal surface, Eq. (33). The probability of 
having a cilium in an area dS around x is denoted by p(x)dS, and we take the 
probability density to be uniform by setting p(x) — 1/.7 for every x of Eq. (33), where 


E d nW? arcsin e 
d 0) sin 0) + r’(0)?dé dp = 14 36 
[ [re sine ior Pa0 de = i+ EAS os 


is the surface area of the spheroid. 
In order to proceed, we now need to balance the drag force with ciliary propulsion, 
and thus need to specify the details of the propulsion mechanism. 


PLOS 


25/73 


@.PLOS | SUBMISSION 


EZ Litostomatea 49 Armophorea 


EC Litostomatea F Colpodea 
HMHM. Spirotrichea EE Nassophorea 


GB Spirotrichea 


Parabasalia 


Ë 
25 E Prostomatea E Armophorea 2 EE Prostomatea E Opalinea 
" EM] Heterotrichea EXT Nassophorea FEES Heterotrichea E Colpodea 
= 20 Opalinea BEES Karyorelictea = HB Oligohymenophorea EE Karyorelictea 
S E Oligohymenophorea S 15 
5 15 5 
= = 
= = 
o ol 
O 10 C 
S 
5 5 
ji M = Bea | = H = Ë 
0 200 400 600 800 1000 1200 0 1500 3000 4500 6000 


Cell length B [um] Swimming speed U [um s7!] 


Figure 19. Histograms of body lengths, B (jum, left), and swimming speeds, U 
(um s^ 1, right), for the 93 ciliates in the database. Ciliates are by far the largest 
organisms in our database, with the average cell length of 

(B) = 194.87 + 207.45 um (n = 91), and an average swimming speed 

(U) = 1147.57 + 1375.64 ums! (n = 81). 
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Figure 20. Histograms of aspect ratios W/ B (left) and body-to-cilium length B/£ 
(right) for ciliates. Most of the cells are prolate, with the mean aspect ratio 

(W/B) = 0.49 + 0.22 (n = 86). The size-to-flagellum length ratios have average values 
(B/€) = 23.13 + 27.03 (n = 26). 


(A) Constant tangential stresses. The simplest model for ciliary propulsion 
assumes that the array of cilia exerts a constant, axisymmetric stress of magnitude + 


along the tangent vector t. Using Eqs. and (35), the total propulsive force can then 


be written as 
Bes | ee emt y, (37) 
S 


with a purely geometric factor given by 
L(B,W) = > | [(r(0) sin 0)? — r(8)r' (0) sin 0 cos 0| dd. (38) 
0 


Balancing the propulsion P, from Eq. with the drag D given by Eq. and 
solving for the swimming speed U; leads then to the theoretical model 


y 
U, = T— ——. 39 
T Seng Bi (39) 


(B) Constant force per cilium. In the second modelling approach, one may 
imagine that each cilium, whose base lies at the point x(0, p), exerts a constant force F 
along the tangent vector t. One cilium then contributes a local thrust along x of 


magnitude 
F 


VORTOL [r(8) sin 8 — r'(0) cos 0]. (40) 


F(—t): e, = 
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If the ciliated cell possesses N such cilia, uniformly distributed over its surface, the 
central limit theorem establishes the total propulsive force to be 


N 
Pr = FL(B, W). (41) 
After balancing with the drag, this leads to the ciliary swimming speed Up predicted by 


this model as 
NT, 


Up = FL——————. 
PF Zorn Dä. SC CFB 


(42) 


(C) Constant surface velocity. The third modelling approach assumes that the 
local speed of the fluid induced by ciliary motion is (almost) constant. To quantify this 
hypothesis, consider a spheroidal cell with a prescribed tangential surface velocity 
distribution us = us(¢)t, where ¢ = cos 0 and t is given by Eq. (34). 

In this case, the Lorentz reciprocal theorem may be used to relate the propulsion 
speed U, of a squirming organism to the surface velocity distribution by 


8 TQ 1 IC 1/2 
u=-Bf (ES) wae (43) 


where 7 = 1/e = 1/4/1 — (W/ B)? is fixed by the morphology of the swimmer. 
Following past work , if we take an almost uniform surface velocity distribution 
of the form 


_ 72 \ 1/2 
Bosak ES | (44) 


where the constant ûs sets the characteristic surface velocity scale, we obtain a model 
for the swimming speed as given by Us = ûs [rë — ro (z — 1) coth `! ro], which may also 
be written in terms of the eccentricity e as 


iu e (1 u =g tanh“! e) . (45) 


With this particular flow assumption, for a very slender cell body (e — 1), Us > tis, 
while for a spherical cell (e — 0) we get U, — 2ti,/3, in agreement with the classical 


result ç 


6.3 Insights from data 


We begin with the constant tangential stress model (A), where the swimming speed is 
given by Eq. (39). In Fig. [21] we plot the measured speed for all ciliated species in our 
database versus the factors accompanying the tangential stress 7 in Eq. (39). The 
scatter of the data points clearly does not support the hypothesis of universal surface 
stress for all organisms. The model can however be used to estimate the effective stress 
T on the surface of each ciliate in the database. The shaded area represents the bounds 
for 7, and fall in the wide range 0.55 — 580 mPa. These values are consistent with the 
estimate of 7 10 mPa for Volvoz colonies [117]. 

A similar comparison for the ‘constant force per cilium’ model (B), quantified by 
Eq. (42), requires the knowledge of the number of cilia N for each swimmer. This 
number is, however, scarcely reported in literature, with only 9 values registered in our 
database. For some species, however, measurements report the number of cilia per unit 
area K, or the distance between neighbouring cilia d. Using the latter, we can estimate 
the number of cilia per unit area to be kq zz 1/d?. Using Eq. (36), K and N can be 
easily related via N = kZ. By doing so, we determined N (equivalently «) for a total 
of 15 ciliated species out of 93, 13 of which had information about the cell swimming 
speed. In Fig. [22] we plot the reported swimming speed versus the right-hand side of 
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Figure 21. The swimming speed U for ciliates plotted versus the numerical factor 
accompanying the constant tangential stress assumed in model (A) and Eq. (39). The 
shaded area encloses all organisms and serves as an estimate of the average effective 
tangential stress for all organisms, with the lower bound of Tmin = 0.55 mPa, and the 
upper bound of Tmax = 580 mPa. Colours distinguish between classes of ciliated 
organisms. The scatter of data suggests that only a large range of values for the stress 
of individual organisms can be inferred. 


Eq. to estimate the effective force per cilium F. We report our estimated values of 
F for each species in Table [2] Our data encloses previous estimates in the range 

0.3 — 1.0 pN [118], and show that the effective tangential forces exerted by each cilium 
may even be two orders of magnitude lower for species like Opalina ranarum. 

The third model (C) assumes the creation of local flows by an almost constant 
surface velocity, whose order of magnitude is fixed by ûs. The predictions of Eq. 
suggest that the swimming speed and the surface velocity are related by a simple 
geometric parameter, namely a function of the cell body eccentricity, e. In Fig. [23] we 
plot the measured ciliate velocities against the theoretical geometric factor determined 
for each species from our data. The model can be used to estimate the magnitude of the 
effective average surface speed for each species. The resulting values span from a few 
tens of um s^! to about 104 ms 1. The average value of the effective surface velocity, 
calculated for all species, (ûs) = 1.42 10? um s^ 1, is about 2 to 3 times the average 
metachronal wave speed we estimate from our data, Amw f, where Àww is the 
wavelength of the metachronal wave created by the collective ciliary beating at 
frequency f. Here also our data confirm and extend previous estimates. For example 
in Ref. [119], tracking microscopy and fluid velocimetry were used to determine with 
precision the flow field of a freely swimming Volvos colony, resulting in estimates of the 


surface speed ûs ~ 200 — 250 uum s^! for species swimming at U ~ 100 — 150 ums. 


7 Conclusion 


7.1 Summary and perspective 


In this paper, based on an initial selection of six seminal papers in the field of biological 
fluid dynamics and physics, we assembled a summary of the experimental data 
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Figure 22. Reported swimming speed U plotted against the numerical prefactor of 
Eq. (22). assuming a constant effective force per cilium in the propulsion model (B). 
Square symbols mark organisms for which the prediction was directly calculated from 
the available data, while circles represent those for which we estimated the number N of 
cilia. Colours distinguish the different taxonomic classes. The visible large scatter of 
data sets the bounds for the effective force per cilium to be in the range of 1.10 107? to 
3.19 pN, represented by the shaded area in grey. 


produced to date on the characterisation of motile behaviour of unicellular 
microswimmers. The material gathered provides a convenient and practical reference 
point for future studies. Our database includes empirical data on the motility of four 
categories of organisms, namely bacteria (and archaea), flagellated eukaryotes, 
spermatozoa and ciliates. Whenever possible, we reported the following biological, 
morphological, kinematic and dynamical parameters: species, geometry and size of the 
organisms, swimming speeds, actuation frequencies, actuation amplitudes, number of 
flagella and properties of the surrounding fluid. In all cases, we also give the 
appropriate references to the publications reporting the measurements. We then 
analysed this information by characterising some of the statistical properties of the cells 
in our database and by introducing theoretical models for each main species in order to 
establish guiding principles for the presentation of the data. We particularly focused on 
the dependence of the swimming speed on the characteristics of the swimmers and 
environmental properties. The analysis shows that qualitative trends established in the 
theoretical framework based on motility in Stokes flows agrees broadly with the 
reported data but that the large degree of variability among species precludes drawing 
general conclusions from the dataset. The modelling approaches can however be helpful 
in rationalising the data, pointing out the relevant dynamic quantities governing the 
locomotion of each individual group. In particular, our data confirm and extend 
estimates of these parameters previously reported in the literature. 

An important result highlighted by our study is that a tremendous statistical 
variability exists in the available data, not only within domains but also within 
smaller taxonomic groups. Little is known about the variability of motility within 
individual species in Nature, neither in terms of their morphological characteristics 
(e.g. size and shape distribution), nor in the details of their propulsion (flagellar or 
ciliary motion). In fact, for every single set in our database, it is not clear at all how 
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Species F [pN] 
Blepharisma sp. 1.25 rz? 
Coleps hirtus 2.89 
Didinium nasutum 8.82 107! 
Opalina ranarum 1.10 107? 
Paramecium caudatum 1.07 107? 
Paramecium multimicronucleatum | 8.04 107! 
Paramecium spp. 2.09 107! 
Spirostomum sp. 2.49 107? 
Stylonichia sp. 2.57 
Tetrahymena pyriformis 3.48 1071 
Uronema marinum 3.19 
Uronema sp. 7.43 107! 
Uronemella spp. 5.78 107! 


Table 2. Estimated values of the effective tangential force F exerted by each cilium for 
the species in Fig. 
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Figure 23. Reported swimming speeds from our database U plotted against the 
geometric factor from Eq. for the constant-flow model (C). The data can be used 
to estimate the range of effective surface velocities to be in the range from 63.0 jum s^! 
to 1.10 10* ms 1, in the grey shaded area. Colours allow to distinguish different 
taxonomic classes. 


representative any particular measurement is of a group of similar organisms in the 
same environment? How sensitive are the propulsion characteristics of these cells to 
changes in the environmental stimuli and how do they adapt to new conditions? With 
the enhanced capacity to process large datasets and with new developments regarding 
automation of image analysis, the task of gathering and processing statistical data is 
becoming increasingly feasible, and new works will be able to discover the fundamental 
principles dictating the locomotion of similar species within the same taxonomic group. 
'The database in its current form, which is stored on the OSF repository (54), would 
benefit from the collaborative effort of the community. By growing further, it would 
help provide up-to-date information on the dynamics of a variety of organisms and 
populations, hopefully further encouraging collaborations between cell biologists and 
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physical scientists. To aid this process, our database is available on GitHub [55], where 
it can be extended and enriched. 

We gave our database the appellation of BOSO-Micro. The first term stands for 
Bank Of Swimming Organisms while the second is there to emphasise that we have 
focused our work on microscopic unicellular organisms. We hope that new versions of 
the database, BOSO-X, will be built by focusing for X on different organisms. An 
obvious suggestion would be to assemble a BOSO-Fish database, given the large amount 
of experimental, computational and theoretical knowledge on the swimming of fish. We 
hope that building exhaustive databases of this sort will further facilitate the work of 
physical scientists on biological problems related to locomotion. 


7.2 Caveats and limitations 


The collection of data gathered in our database is inevitably incomplete and biased, in 
particular due to the way the initial set of literature sources, focused on biophysical 
studies, was chosen. Despite our efforts to carry out a broad search for swimming data, 
it is possible that important references were left out. The mitigation strategy in this 
case relies on making the database public and expandable [55]. 

Regarding the presented data, a major limitation is of course their sparseness. The 
relevant parameters in the description of motility are incomplete for many species, 
especially the variables related to the beating of cilia and flagella, which hinders direct 
comparison with theoretical models. 

Furthermore, the database was populated using data presented across different 
papers, books, registers and reports, and the multiplicity of sources introduces a 
significant and inherent noise. For many species, reported measurements of one or more 
characteristics refer to different experimental environments. Even if those are reported, 
different strains of the same species may behave differently under slightly modified 
physical and chemical conditions. 

It is also important to highlight the limitations and assumptions of the models used 
in our paper. The models were designed to assist the presentation of data in the context 
of established ideas regarding microscale locomotion, and to provide quick estimates of 
the relevant dynamic characteristics of microswimmers. Common to all the models is 
our assumption that the cell bodies are spheroidal. A look at Fig. [3] quickly reveals that 
this hypothesis is a crude approximation for many species in our analyses (e.g. 
Caulobacter crescentus, Ceratium tripos, Stentor). We have made this choice in 
modelling in order to account for the influence of both the cell body length and width 
in an analytical way. The diversity of form, which might be crucial for certain 
locomotion strategies, has no reflection in the considered simplistic models, yet it must 
be incorporated into specific models describing particular organisms. Similarly, in the 
case of swimming eukaryotic cells, several of our hypotheses on the flagellar beat ought 
to be examined carefully. For spermatozoa and flagellated eukaryotes, we assumed the 
form of a simple sinusoidal wave, whereas many species display more complex flagellar 
beating patterns (e.g. complex waves displayed in Columba livia and Sturnus vulgaris 
spermatozoa). For flagellated eukaryotes, we have neglected hydrodynamic interactions 
between flagella, which is a simplified approximation. In the case of ciliates, the three 
models we have introduced also do not take into account hydrodynamic interactions 
between neighbouring cilia, nor the effect of the polarised beating of cilia and their 
recovery stroke. Despite these limitations, we hope that the use of modelling may also 
prove useful in rationalising and organising future data on swimming organisms along 
similar lines. 
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A The database of swimming microorganisms 


In Table [3] we present a short glossary with the main symbols used in the database. 


Symbol Meaning Unit 
B Body length pm 
W Body width um 
N Number of flagella or cilia - 

L Lengths (mostly flagella, otherwise specified) um 
Nw Number of waves (full periods, or crests) produced by flag- - 
ellar beat 
À Wavelength of flagellar waves (of helicoidal body and of pm 
metachronal waves indicated by a subscript B and MW, 
respectively) 
A Length of a complete wave along the flagellum (or path, pm 
indicated by subscript) 
h Amplitude of waves (for helicoidal bodies, a subscript B um 
added) 
U Swimming velocity ums 
Ww Flagellar beat frequency a) 
Q Frequency of the rotation of cell body s71 
c Wave speed of flagellar beat (or metachronal wave) ums! 
V Volume of cell body um? 
£ Length of cilia pm 
d Distance between cilia pm 
b Radius of flagella um 
K Number of cilia per unit area um ? 
f Beating frequency of cilia s71 
G Gyration (frequency at which organisms revolve around the — 


axis of movement) 
"n Viscosity of the swimming medium mPas 
Table 3. List of symbols used in the database, together with their explanation and 
units. 


For every entry in the database, in the case when more than one measurement was 
available, we report the average value and the standard deviation using the + notation. 
Values inside parentheses specify the range of the values measured, e.g. (min — Tmax). 
Sometimes only the upper boundary was available, indicated by a preceding ‘max’. 
When the information was not available in the texts of the articles, the figures or the 
graphics were analysed with the GNU Image Manipulation Program (GIMP) software 
in order to extract data. This is indicated in the tables by a superscript * or * 
respectively, if figures or graphics were used. 

The various tables of data are organised as follows. Table [4] contains the data for 78 
organisms in the branch of bacteria (with 5 spiral-shaped bacteria included). 
Spirochaetes (18 species) and Spiroplasma (2 species) were separated from the other 
bacteria because of their distinct mode of locomotion and are presented in Table [5] The 
data for the 10 species of archaea are contained in Table [6] 

Eukaryotes have also been divided into three groups. The data for flagellated 
eukaryotes (121 species) are presented first in Table [7] followed by spermatozoa (60 
species) in Table [B| and finally ciliates (93 species) in Table [9] 
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